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We predict that pairs of polar molecules can be weakly 
bound together in an ultracold environment, provided that a 
dc electric field is present. The field that links the molecules 
together also strongly influences the basic properties of the re- 
sulting dimer, such as its binding energy and predissociation 
lifetime. Because of their long-range character these dimers 
will be useful in disentangling cold collision dynamics of po- 
lar molecules. As an example, we estimate the microwave 
photoassociation yield for OH-OH cold collisions. 



A recurring theme that has emerged in the study of ul- 
tracold matter has been one of control over interparticle 
interactions. Most notably, magnetic fields and res- 
onant laser light || have been employed to alter atomic 
scattering lengths, influencing in turn the behavior of 
Bose-Einstein condensates. Perhaps even more intrigu- 
ing is the ability, using the same tools, to catalyze the 
production of diatomic molecules [rj-pdf| . When the ini- 
tial atomic sample is Bose condensed, there can result 
coherent superpositions of atoms and molecules [jTl| in a 
process dubbed "superchemistry" Q . More broadly, ex- 
ternal fields may be expected to guide chemical reactions. 
Chemical processes involving atom exchange have been 
considered at ultralow temperatures |L3|-|lq|, and have 
been observed so far experimentally in the context of 
three-body recombination p],pJ| and vibrational quench- 
ing§. 

A prerequisite to effective control over atoms or 
molecules is a clear understanding of their interactions at 
ultralow temperatures. Even for alkali atoms, ab initio 
potentials have insufficient accuracy to predict scattering 
lengths; this statement is even more true for molecules. 
Investigators must therefore rely on key experimental 
data that will unravel threshold scattering properties of 
molecules in the most straightforward way. For the alkali 
atoms, an efficient tool for extracting scattering lengths 
has been optical photoassociation (PA) to weakly-bound 
states of alkali dimers. The states in question are gener- 
ated by spin-orbit couplings that lead to avoided cross- 
ings in potential curves at large values R of interatomic 
separation (typically tens of Bohr radii) raj. Because of 
their long-range character, they are insensitive to details 
of the small- R potentials that are poorly characterized. 
Free-to-bound optical spectra of these states therefore 
contain, in a reasonably accessible form, important in- 
formation on scattering phase shifts |18|. 

More recently, even more exotic molecular states 



have been predicted to exist in ultracold environments. 
Greene and collaborators have predicted Rydberg states 
of diatomic molecules with elaborate electronic wave 
functions, whose density plots resemble trilobites [51| 
or butterflies p5| . In addition, Cote, Kharchenko, and 
Lukin have identified possible molecules composed of 
hundreds of neutral atoms clinging to a single impurity 
ion in a Bose-Einstein condensate fl2jj . 

In this Letter we introduce another set of unusual 
molecular states. These are composed of two ground- 
state polar molecules held at large intermolecular sepa- 
ration under the joint influence of electric dipole forces 
and an external electric field. The properties of these 
metastable dimers, such as binding energy, intermolec- 
ular separation, and predissociation lifetime, depend 
strongly on the value of the field, leading to new types 
of control over intermolecular interactions. These states 
may be populated by rapid electric field ramps, or else 
by the analog of photoassociation spectroscopy in the mi- 
crowave regime. Below we describe the origin and proper- 
ties of these states, and estimate the trap loss rates asso- 
ciated with their formation and subsequent dissociation 
in microwave PA experiments. In this way we sketch a 
general strategy for understanding and manipulating cold 
collisions of polar molecules produced by either Stark de- 
celeration |22| or buffer-gas cooling |23| . 

For concreteness we will consider dimers of two OH 
molecules, to be denoted [OH]^, although the phe- 
nomenon should apply quite generally to polar molecules, 
for example NH3 |^2| , that suffer a linear Stark effect at 
experimentally realizable fields. We have previously de- 
scribed a model for the OH-OH interaction Hamiltonian, 
including the parity, spin, and nuclear spin degrees of 
freedom that are most relevant at ultralow temperatures 
|^| . The OH ground state has spin j = 3/2, which, cou- 
pled to the nuclear spin i = 1/2 of the proton, gives total 
spin / = 2 in the molecular states we consider. For clar- 
ity, we restrict the present model to molecules that are 
spin polarized with rrif = f = 2, and neglect couplings to 
different values of mj. Further, we incorporate only the 
partial waves L = 0, 2. These restrictions do not change 
the qualitative behavior of the model p4j, but they do 
render its structure more tractable. 

The scattering thresholds in this model are set by the 
Stark splitting A(£) of the molecular levels. At the rela- 
tively low fields of interest here, rotational states are not 
substantially mixed by the field, but the opposite parity 
states of the molecule's Lambda doublet are mixed. Thus 
there is a weak tendency for the field to align the dipole 
moments, along the field when uj < 0, and against the 
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field when lo > 0, where u> is the projection of j along the 
body axis of the molecule (these alignments are reversed 
when the laboratory projection rrij < 0) p5|] . For con- 
venience we refer to these alignments, valid at infinitely 
large intermolecular separation, as f and j, respectively. 
A set of adiabatic potential energy curves for our model is 
shown in Figure 1(a), where the dc electric field strength 
is taken to be £ — 10 4 V/cm. 
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FIG. 1. (a) Adiabatic potential energy curves for the 
six-channel model of [OH]2, as described in the text. Avoided 
crossings in these curves create the potential well labeled 
Vb, which holds bound states of the dimer. (b) shows the 
electric-dipole coupling matrix element between Vb and the 
ground state Vq. 

When the separation between the molecules is not in- 
finite, the partial wave L is no longer a good quan- 
tum number, different values being mixed by strong 
anisotropies of the dipole-dipole interaction. The adi- 
abatic channels shown therefore represent different ori- 
entations of the intermolecular axis relative to the 
laboratory- fixed electric field axis, as indicated schemat- 
ically by the diagrams within ellipses. For example, con- 
sider the highest threshold where both molecules point 
against the field. The attractive adiabatic potential cor- 
responds to a combination of L = and 2 partial waves 
that places the dipoles end-to-end on average, whereas 
the repulsive curve places the dipoles side-by-side on av- 



erage. 

When the molecules approach still closer to one an- 
other, there is an inevitable curve crossing between the 
repulsive potential correlating to the | f }+ | J.) threshold 
and the attractive potential correlating to the | |)+ | |) 
threshold. These channels are dipole-coupled, leading to 
an avoided curve crossing at R = 60 a.u. in the figure. 
The resulting potential curve, labeled VJ, in Figure 1(a), 
is deep enough to hold a set of quasi-bound states of the 
[OHjj dimer. (Note that this curve was denoted "C/ " 
in Ref. [Q. Here we label it more consistently with the 
theory of photoassociation, to which we refer below.) 

The potential V&, and its bound states, are strongly 
subject to the strength of the applied electric field £ . 
First, changes in £ shift the relative thresholds, chang- 
ing the location of the crossing point. Second, the cou- 
pling matrix element between the channels that cross 
to form Vb varies with field as ~ k/(l + k 2 ), where 
k = 2/i£ /A(0), and /i and A(0) are the molecule's dipole 
moment and zero-field Lambda-doublet splitting, respec- 
tively ^4j. This means that in low field the crossing is 
purely diabatic, to lowest order, and the molecules can 
approach to small values of R where hydrogen bonding 
and chemical forces can act. 

That the field controls these states is illustrated in Fig- 
ure 2. Here the binding energies of potential Vb are plot- 
ted as a function of electric field £, from zero to 10 4 
V/cm. The solid lines refer to a close-coupling calcula- 
tion in the six-channel model, whereas the dashed lines 
give the bound state energies of the adiabatic potential 
Vb. Similar results are obtained in the more complete 
calculation described in Ref. [24|| . A striking feature of 
this figure is that for fields below ~ 1000 V/cm there are 
no bound states at all, because the potential becomes too 
shallow. Because of the essential role the field plays in 
creating these states, we refer to them as "field linked" 
states. 

The field-linked (FL) states remain coupled to chan- 
nels with lower-energy thresholds, and can therefore pre- 
dissociate to these channels. For several selected field 
values the energy widths for predissociation, computed 
from the eigenphase sums in the close- coupling calcula- 
tion, are indicated by error bars on the energies. Two 
general trends are readily apparent: first, the widths as a 
whole become exponentially narrower as the binding en- 
ergy approaches threshold; second, there is an additional 
field-dependent modulation in the widths, governed by 
the varying overlap integrals between the resonant states 
and the continua into which they decay. 

The dominance of the electric field over the FL states 
implies that ultracold molecular collisions can be ma- 
nipulated. For example, resonant control over scattering 
lengths in | |}+ | |) collisions is possible by tuning bound 
states of Vb through the scattering threshold [Q . A sim- 
ilar tuning was noted for alkali atoms in extremely strong 
electric fields [27[] , but without making the explicit con- 
nection to bound states. 

The FL states also imply a degree of control over chem- 
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ical reactions. For example, placing two OH molecules in 
a FL state holds them at a "safe" distance where the 
reaction cannot occur (assuming for example that the 
exoergic reaction OH + OH — » H2O + O actually occurs 
at ultralow temperatures). The reaction can then be ac- 
tivated at will by switching off the field, which could 
be done at any desired rate simply by varying the field 
strength. 
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FIG. 2. Binding energies of bound states in the potential 
Vb versus electric field. Solid lines: close-coupling results, 
dashed lines: single-channel adiabatic appoximation. Error 
bars: approximate predissociation widths, as computed form 
the eigenphase sums in the close-coupling calculations. 

FL states can presumably be created from a gas of | [) 
molecules by quickly ramping the electric field, just as 
bound states of Rb2 have been created in magnetic field 
ramps . Setting the field to an appropriate value will 
guarantee that the state created will be long-lived against 
predissociation. It is also possible that they have been 
produced already by collisions in the supersonic expan- 
sion that feeds the Stark slowing process |22|] , although 
they have not yet been sought in this environment. 

As a practical issue, a gas of molecules in the weak- 
electric-field-seeking | I) state is subject to large colli- 
sional losses in an electric field and is therefore not suit- 
able for evaporative cooling J24|. The | f) state, on the 
other hand, is collisionally stable, being the lowest energy 
state of the OH molecule, and it can be trapped in either 
a magnetostatic or an off-resonant optical dipole trap. 
Adiabatic curves correlating to the | f)+ I T) threshold 
do not experience a curve crossing at large i?, and allow 
the molecules to approach to small values of R (see Fig. 
1). The cold collision parameters, for instance the scat- 
tering length, will then have to be measured in a suitable 
experiment. 

We are therefore led to consider photoassociation (PA) 
spectroscopy of the molecular pair. This will be PA at 
microwave, rather than optical, wavelengths, appropriate 



to the energies of the FL states. We follow here the dis- 
cussion of Ref . |2(| , and consider single-channel processes 
only. This will enable order-of magnitude estimates of the 
multichannel response. The PA process may be viewed 
as a scattering event, where two molecules incident in 
channel Vq with relative energy E are promoted by a mi- 
crowave photon of detuning A to a FL state in potential 
Vb- This state can then predissociate, typically releas- 
ing a large fraction of the photon's energy and leading 
to trap loss in the exit channel. The rate constant for 
collisional trap loss via a bound state in potential V& is 
given by [§6) 



«■ k 2 (£_ A )2 + (( 7 + r)/2) 2 ' v ' 

where v and k are the relative incident velocity and wave 
number, 7 is the predissociation width, and the stimu- 
lated width r is given by the Fermi's Golden Rule ex- 
pression 

r = 27r|(/ o |y rad |0 fc }| 2 . (2) 

Here K-ad — < ~M'£ac > is the independent dipole matrix 
element for the transition, averaged over all coordinates 
except R |28| . £ ac is the amplitude of the ac electric field 
of the applied microwaves, which we assume are polarized 
parallel to the dc field. The width ([!]) also depends on 
fo(R), the energy- normalized incident wave function, and 
4>b(R), the unit-normalized bound state wave function in 
potential Vb- On resonance, and in the typical case r <C 7 
the rate is approximately 

Anv r 

k 7 

hence is proportional to the microwave intensity. 

Measurements of optical PA spectra have proven use- 
ful in assessing the scattering lengths of alkali atoms [^9| , 
largely because the overlap integral in (||) is sensitive to 
the amplitude of the wave function /o at the intermolec- 
ular distance R c where the transition occurs. In fact, 
this circumstance enables PA rates to be estimated using 
a "reflection approximation" that focuses on the wave 
functions at R = R c alone |3Cj . The FL states should ex- 
hibit a similar sensitivity to scattering length, although 
now the independence of the dipole coupling V la d cannot 
be neglected, as seen in Figure 1(b). In the limit of large 
R, ^rad tends to zero, owing to the inability of a single 
photon to drive a transition of both molecules at once. 
At smaller values of R, where the avoided crossing occurs 
that creates Vb, the radiative coupling attains a nonzero 
value. 

To illustrate the feasibility of PA experiments in this 
system, and to emphasize the strong sensitivity of the 
spectra to electric field and scattering length, we estimate 
rate constants for several different circumstances. These 
results are summarized in Table I. In all cases we have 
assumed that the ac microwave field has a field strength 



3 



of £ ac = 15 V / cm, which can be reasonably achieved in 
the laboratory . We have computed the resonant PA 
rate according to Eqn. ([!]) , assuming two different short- 
range potentials that produce the scattering lengths a 
til) a.u. and a — —400 a.u. in the incident Vq ground 
channel, in zero electric field. The collision energy is 100 
fiK in the results shown here. 

The Tables show, for several values of the dc elec- 
tric field £, the binding energy of several states; their 
predissociation widths 7; their stimulated widths T for 
each of the scattering lengths; and the resonant loss rate 
K\ oss , again for both scattering lengths. The table illus- 
trates that a fairly small field change, on the order of 
1 kV/cm, can completely alter the loss rates. Moreover, 
the rates are different for different scattering lengths, im- 
plying sensitivity to details of the short-range part of Vq. 
A complete map of loss versus field ought therefore to 
contain a wealth of information about cold collisions of 
OH molecules. Equally importantly, the absolute mag- 
nitudes of the loss rates can be quite large, as high as 
several 10 -11 cm 3 /sec, which should be easily measur- 
able if densities are high enough for molecular collisions 
to be observed at all. 

In summary, we have introduced a set of weakly-bound 
"field-linked" states of ultracold polar molecules that ex- 
ist only in the presence of an electric field. The sensitivity 
of these states to electric field suggests an unprecedented 
level of control over intermolecular forces, since the very 
shape of the interaction potential can be manipulated. In 
addition, their purely long-range character makes them 
invaluable tools for approaching the difficult problem of 
deciphering ultracold molecular collisions. 

This work was supported by the NSF. We acknowledge 
useful discussions with J. Hutson and J. Ye. 
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TABLE I. Characteristics of field-linked states and the 
photoassociation yield subject to an £ = 15V7 cm ac electric 
field, at a collision energy E = 100 fiK. Numbers in paren- 
theses denote the power of 10 by which the numbers are to 
be multipled. 
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